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This is the first of three papers devoted to the nonequilibrium thermodynamics of amorphous 
materials. Our focus here is on the role of internal degrees of freedom in determining the dynamics 
of such systems. For illustrative purposes, we study a solid whose internal degrees of freedom are 
vacancies that govern irreversible volume changes. Using this model, we compare a thermodynamic 
theory based on the Clausius-Duhem inequality to a statistical analysis based directly on the law of 
increase of entropy. The statistical theory is used first to derive the the Clausius-Duhem inequality. 
We then use the theory to go beyond those results and obtain detailed equations of motion, including 
a rate factor that is enhanced by deformation-induced noisy fluctuations. The statistical analysis 
points to the need for understanding how both energy and entropy are shared by the vacancies and 
their environments. 
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I. INTRODUCTION 

This is the first of three papers describing our ef- 
forts to develop a thermodynamically well-founded the- 
ory of nonequilibrium phenomena in amorphous materi- 
als. Specific goals of this project are to develop a thermo- 
dynamic understanding of the effective disorder temper- 
ature and the role that it plays in Shear- Transformation- 
Zone (STZ) theories of amorphous plasticity [U, [2, S S Hi ■ 
While working toward these goals, we have encountered 
a number of fundamental questions. Those questions 
include: What is the most basic statement of the sec- 
ond law of thermodynamics? How can we reconcile the 
different approaches to nonequilibrium thermodynamics 
taken by engineers, applied mathematicians, and physi- 
cists? Are the dynamic roles played by internal degrees 
of freedom properly described by any of those theoretical 
approaches? Many recent developments in the physics of 
glassy materials, including the STZ theory, are based on 
the idea that the state of disorder in such systems is de- 
scribed by an effective temperature that is not necessarily 
the same as the ordinary temperature. In what sense is 
the effective temperature a well defined thermodynamic 
concept? 

In this first paper, we focus on questions regarding 
internal degrees of freedom. For illustrative purposes, 
we address that issue in the limited context of a simple 
model of a uniform, not necessarily glassy, solid in which 
vacancies govern irreversible volume changes. We origi- 
nally developed this vacancy model as a way of studying 
irreversible changes in the volume of a glassy material 
subject to varying temperatures and pressures. We hope 
to return to such applications in the future; but, for the 
present, we use the model purely as an aid for explor- 
ing theoretical ideas. In the second paper [Q], we use the 
insights gained here to define the effective disorder tem- 
perature and to write equations of motion for it. Finally, 
in the third paper 0] , we reformulate the STZ theory in a 



way that is consistent with the thermodynamic analysis 
presented here and in Q. 

Almost by definition, the irreversible responses of ma- 
terials to applied forces are determined by internal de- 
grees of freedom. Deforming amorphous solids are gen- 
erally described in terms of internal entities such as flow 
defects or, of special interest here, STZ's. Qualitatively 
similar internal structures appear in theories of disloca- 
tion motion in crystalline solids, and in nonequilibrium 
theories of granular materials and complex fluids. The- 
ories of these dissipative phenomena necessarily invoke 
the second law of thermodynamics, at the very least as 
a constraint on the equations of motion for the internal 
variables. 

There is a very large body of literature on this sub- 
ject. For example, see monographs by Lubliner [8j], Mau- 
gin [9], and Nemat-Nasser which we have found to 
be especially useful. Essentially all of this literature is 
based on the postulate that the Clausius-Duhem entropy- 
production inequality is the fundamental statement of 
the second law of thermodynamics; therefore, we refer 
to that theoretical starting point as "conventional" . We 
recognize, however, that the body of literature to which 
we are referring contains many different points of view, 
and that these points of view have continued to evolve in 
recent decades, especially in the engineering and applied 
mathematics communities. Our motivation for develop- 
ing a statistical approach based directly on the law of 
increase of entropy simply reflects the fact that we have 
not been able to take any version of the conventional ap- 
proach far enough to answer the questions that we are 
asking. 

In this paper, we look at the issues concerning the sec- 
ond law and internal degrees of freedom in a simple but 
physically realistic situation ~ the vacancy model men- 
tioned earlier. We start in Sec. |lT] by introducing the 
model and then, in Sec. IIIIl by briefly summarizing a 
conventional analysis. The first-principles statistical the- 
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ory and the resulting nonequilibrium equations of motion 
are presented in Sees. |IV]-|Vll We conclude in Sec. IVIII 
with some remarks about the broader implications of our 
results. 



II. VACANCY MODEL 

We consider a uniform viscoelastic solid of total volume 
V, containing a small but extensive number of vacancies 
Ny. It may be easiest to visualize this solid as being 
noncrystalline, but that assumption is not essential for 
present purposes. To avoid the complications of position 
dependent deformation, we assume that this system re- 
mains spatially uniform at all times, and we work with 
extensive quantities rather than local densities. We fur- 
ther assume that the system is never too far from thermo- 
dynamic equilibrium, i.e. that we are not dealing with ex- 
tremely rapid nonequilibrium phenomena for which local 
thermodynamic concepts would be invalid. This quasi- 
equilibrium condition is essential for our arguments in 
Sec. lYl 

Our model, like any model of a material subject only 
to volume and not shape deformations, is technically vis- 
coelastic rather than elastoplastic. If it is allowed to equi- 
librate at a fixed temperature and pressure, it eventually 
returns to the same equilibrium volume with the same 
number of vacancies; whereas a true elastoplastic ma- 
terial would return to a permanently deformed shape if 
subjected to shear. A second difference is that a plastic 
material can undergo steady-state shear flow, and usu- 
ally exhibits a yield stress that marks the onset of that 
behavior. No such steady-state deformation can occur in 
the present case. 

The fundamental differences between these kinds of in- 
elastic deformation are important, but are not the cen- 
tral issues to be discussed here. Rather, the model of 
purely volumetric deformation is especially useful to us 
because there is no need for anything other than an addi- 
tive decomposition of arbitrarily large elastic and inelas- 
tic volume changes. No special mathematical efforts are 
needed to compute the results of complex sequences of 
deformations; but the separation between elastic and in- 
elastic deformations remains a nontrivial topic of interest 
as discussed, for example, in plj . 

The volume V in this model consists of three additive 
components: 

V = Vo + Vel+V,n. (2.1) 

Here, Vq is a reference volume, determined by the entropy 
(or temperature). For simplicity, we neglect thermoelas- 
tic effects and assume that Vq is just a constant. The 
elastic volume, Vei, is associated with reversible changes 
in the elastic energy. An increment SVei is a change in 
the total volume at fixed entropy and fixed iV„; that is, 
it takes place with no change in the internal state of the 
system. Our central assumption is that the inelastic vol- 
ume associated with the vacancies is simply Vm = vq Ny, 



where vq is the effective volume of a vacancy. Vei and 
Vin are independently "variable" but not independently 
"controllable." In equilibrium, Vei is controlled directly 
by the pressure. On the other hand, Ny is a "hidden" 
internal variable so long as the system is not coupled to 
a chemical-potential reservoir that controls the number 
of vacancies. Nevertheless, we must treat Ny as having 
its own dynamics, and being able at any time to change 
in ways that are not directly constrained by concurrent 
changes in Vei or the entropy. Specifically, we assume 
that the rate of inelastic volume deformation is 

Vn^voNy, (2.2) 

and that Ny is a dynamical variable that satisfies its own 
equation of motion. 

III. CONVENTIONAL THEORY 

A conventional analysis of this model, described along 
lines laid out by Coleman, Noll and Gurtin in the 1960's 
starts by writing the first law of thermodynamics 
in the form 

tl^-pV + Q, (3.1) 

where U is the internal energy, p is the pressure, —p V is 
the work done on the system and Q is the rate at which 
thermal energy is entering the system. The conventional 
theory then postulates that there exists an entropy S and 
a temperature 9 defined by a continuity equation. For 
this spatially uniform system, that equation is simply 

^-f-S. (3.2) 

Here, the temperature 9 is expressed in energy units 
(fcs — 1); Q/9 is the rate at which entropy is entering the 
system; and S is the entropy source, i.e. the rate at which 
entropy is being produced. The conventional statement 
of the second law, the Clausius-Duhem inequality, says 
that the entropy production rate is non-negative 

E > 0. (3.3) 

These relations are taken to be axiomatic; they do not 
presuppose any statistical interpretation of S or 9. 
Eliminating Q between Eqs. p.ip and (|3.2p . we find 

9S -if -pV = 9J:>0, (3.4) 

which is conveniently rewritten by transforming to the 
Helmholtz free energy F(6l, Vei, 7V„) = U{S,Veu Ny)~9 S , 
and then performing the partial differentiations 

= -01] < 0. (3.5) 
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Here we have used Eqs. (12.11) and (|2.2I1 . 

In the spirit of Coleman and Noll ^iJl, we recognize 
that this expression consists of three separate, indepen- 
dent inequalities because, as discussed above, the time 
derivatives are unconstrained by each other. We satisfy 
the first inequality by identifying S = —dF/dO, thus 
recovering the familiar thermodynamic relation. The 
second inequality usually is satisfied by identifying p — 
— dF/dVei i.e. using the equilibrium relation between p 
and Vei- However, if we were interested in thermovis- 
coelastic effects, then we would satisfy this inequality by 
writing a dissipative equation of motion for Vei of the 
form 



Vel = -7e 



F,i ^F+pVel 



(3.6) 



where jei is a non- negative dissipation coefficient. In 
this way, we would account for the energy dissipation 
that accompanies the relaxation of the "viscous pressure" 
p + dF/ dVei ■ Since thermoviscoelasticity is not the topic 
of primary interest here, we simply adopt the equilibrium 
relation from here on. 

Finally, using Eq. (|2.2p . we write the third inequality 
in the form 



dF, 



Ny>0, Fy=F+pvQNy. (3.7) 



This is a specific realization of the Clausius-Duhem in- 
equality, Eq. (|3.3p . The same result was obtained in 
Rice's classic 1971 paper [l^l where, however, the free 
energy F was assumed to be a function of the total de- 
formation rather than the elastic part alone, so that the 
pressure was missing in the expression for Fy. The in- 
equality of p.7p is satisfied by 



dF, 
dN„ 



(3.8) 



where "fy is again a non-negative dissipation coefficient. 
More generally, any monotonically increasing function of 
Ny that vanishes where dFy/dNy — can be used on the 
right-hand side of Eq. ([3^ . 

Ever since Coleman and Noll introduced their ax- 
iomatic version of thermomechanics, physicists (such as 
JSL) have been impressed by its mathematical elegance, 
but have worried that it might be incomplete because it 
does not start with a statistical definition of entropy. It 
is not clear what statistical interpretation of the entropy 
S is implied by the preceding equations or, conversely, 
how the internal energy U might depend on S. For ex- 
ample, it is not obvious in a conventional formulation 
how to evaluate the free energy Fy in Eq. p.7p . More 
importantly, the Coleman-NoU postulates operationally 
define a temperature as well as an entropy. In the anal- 
ysis presented here, we are looking ahead to an effective 
temperature theory in which there will be two different 



temperatures - a situation that seems to be beyond the 
scope of the conventional axiomatic formulation. 

Lastly, we note that the axiomatic approach makes no 
mention of a thermal reservoir. It seems to us that any 
theory of this kind ought to include a specific mechanism 
by which the temperature is controlled. If that mecha- 
nism involves coupling to a thermal reservoir, then the 
theory ought to predict the rate at which heat is flowing 
between the system and the reservoir. Conversely, the 
theory should be able to predict what happens if that 
flow is constrained, as in an adiabatic process. But any 
coupling to a reservoir disappears when Q is eliminated 
in Eq. 



IV. STATISTICAL THEORY 

The basic statistical statement of the second law is 
that the system as a whole, including any thermal reser- 
voir to which the subsystem of primary interest may be 
coupled, must move toward states of higher probability, 
i.e. to states of higher entropy. Although the Coleman- 
NoU procedure assigns no a priori statistical significance 
to the entropy, this principle lies at its heart. In their 
formulation, however, the principal focus is on spatial 
heterogeneities. The entropy of the system as a whole 
increases as heat flows between spatially separated ele- 
ments, each of which is always in a state of local equi- 
librium with its own local energy, entropy, and tempera- 
ture. It is conceptually easy, albeit mathematically more 
complicated, to add spatial heterogeneity to the vacancy 
model. We do not do this explicitly here but, neverthe- 
less, anticipate the need to reinterpret our uniform model 
as just one element of a larger, spatially inhomogeneous, 
coarse-grained system. 

Our strategy is to start with a statistical deflnition of 
entropy, and to introduce a thermal reservoir, but oth- 
erwise to stay as close as possible to the conventional 
analysis. Therefore, in analogy to Eq. p.ip . we begin by 
writing the first law in the form 



pV^U + Ur, 



(4.1) 



where Ur = Ur{Sr) is the energy of the reservoir as a 
function of its entropy Sr. Similarly, in analogy to Eqs. 
p.2p and (|3.3p . the second law is 



Sr > 0, 



(4.2) 



where Sneq is the entropy of a system that is not neces- 
sarily in thermal equilibrium. 

The main question is what to use for Sng^q. We pro- 
pose, with several conditions to be listed below, that the 
correct choice of this entropy has the form 

Sn^giU, V, {A„}) = In niU, v., {A„}), (4.3) 

where i^{U, V, {A„}) is a constrained measure of the num- 
ber of states of the system with energy [/, volume V, and 
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specified values of a set of internal variables {Aq,}. The 
Ac's are out of equilibrium if their values are not the 



ones that maximize Sn 



When all of them do maxi- 



mize Sneq, i-G. when Aq, = A'^ , then we require that 
the equihbrium entropy Seq{U, V) — In il{U, V) be accu- 
rately approximated by 

,(c/,y)«^^„e,(c/,v^,{An). (4.4) 



ySeqy 



This approximation must become an equality in the ther- 



modynamic limit, V 



oo. In general, Sneq < Seq, be- 



cause the constrained entropy Sneq counts fewer states 
than the unconstrained entropy Seq ■ We require that the 
difference between these quantities per unit volume be- 
come negligibly small as {Aq}— >{Aq'} and as the size of 
the system becomes indefinitely large. Without this con- 
dition, we would not have a single, well defined entropy 
upon which to base a self-consistent thcrmomechanical 
theory. 

Validity of Eq. (|4.4p therefore requires that three con- 
ditions be satisfied: 

(1) The set of variables {A^} must be sub-extensive. 
If there are Na such variables, and there are N total de- 
grees of freedom in the system, then Na/N must vanish 
in the thermodynamic limit. Otherwise, the variations 
of the Aq's would produce an extensive entropic correc- 
tion to the equilibrium free energy, and Eq. (|4.4p would 
not be correct. More explicitly, note that we can com- 
pute Seq by integrating over each of the variables Aq 
in fl{U, V, {Aq}), obtaining a correction to In propor- 
tional to TVq. That correction must be negligibly small 
compared to Sneq{U, V, {Aq'^}), which is proportional to 
N. 

(2) We must be working in the quasi-equilibrium limit, 
where all the unconstrained degrees of freedom have 
rapidly come to equilibrium, and where their fluctuations 
have been accounted for in computing Sneq- 

(3) Condition (1) requires that the Aq's be coarse- 
grained variables. If there is only a sub-extensive number 
of these variables, then each of them must be a sum over 
a statistically large number of degrees of freedom. That 
is, the Aq themselves must be extensive. (Of course, 
nothing prevents us from interpreting them as spatial 
averages over a macroscopically large system.) It then 
follows that the entropies associated with each of the Aq 
must be included explicitly in Sneq- For example, our sin- 
gle internal variable Ny describes an extensive population 
of vacancies. The associated entropy, i.e. the logarithm 
of the number of ways in which the Ny vacancies can be 
arranged in the volume V, must be contained in Sneq- 

Accordingly, the entropy appearing in Eq. (|4.2p is 

SneqiU, Vel,Ny) - In Q{U, Vel,Ny). (4.5) 

For reasons discussed in Sec. [TTl we replace V by Vei 
as an independent argument of Sneq- We then invert 
SneqiU, Vel , Ny) to obtain U{Sneq, Vei,Ny). We identify 



and, as stated following Eq. (|3.6p . we use the equilibrium 
thermodynamic relation for the pressure 



dU 
dVei 



(4.7) 



Therefore, 



dU 



U^-pVel + 

The first law, Eq. (|4.ip . becomes 
/ dU 



Ny+eSneq- (4.. 



pVtr 



\dNy 



Ny-UR = 9Sneq, (4.9) 



where we have used V = Vei + Vin to eliminate Vei - 

At this point, we depart from the strategy that led 
to Eq. (13. 4|) . Instead of eliminating the coupling to the 
thermal reservoir as was done there, we use Eq. (j4.9p to 
evaluate Sneq in the second law, Eq. (|4.2p : and we iden- 
tify Sr — Ur/Or, where 9r = dlJR/dSR is the reservoir 
temperature. We also use Eq. (|2.2p to eliminate Vin in 
favor of Ny. The result is 



W{p,Ny,Ny) 



1 



^ I C/fl > 0, 



where 



V^{p,Ny,Ny) 



pvo 



f dU 



XdNy 



(4.10) 



Ny (4.11) 



is the rate at which inelastic work is done on the system 
minus the rate at which energy is stored by the vacancies. 
As will be seen, W is a dissipation rate that appears in 
various forms throughout this series of papers. 

The appearance of Ur in this inequality is important, 
because we control the temperature of the system by con- 
trolling the reservoir temperature. Thus the inequality 
in Eq. (|4.10p must be satisfied for arbitrary variations 
oIUr, independent of whatever else is happening in the 
system. We also must satisfy this inequality for arbitrary 
variations of Ny. For example, the vacancy population 
could be relaxing toward an equilibrium value while Ur 
remains constant. Therefore, in the spirit of Coleman 
and Noll, we argue that the only way to satisfy this com- 
bined inequality for all possible variations of the system 
is to enforce two separate, independent inequalities: 



and 



W{p,Ny,Ny)>0, 



1--^]Ur>Q. 



(4.12) 



(4.13) 



The first of these relations is essentially identical to the 
Clausius-Duhem inequality in Eq. p.7p . The differences 
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are that we have derived Eq. (|4.12p from statistical first 
principles rather than postulated it, and that we know 
exactly what energy and entropy are involved in it. 

The second inequality is satisfied by requiring that Ur 
be a function of 9 that changes sign only when 9 = On', 
therefore we write 



-Un = Ai9,9R)i9n-9) = ^ 



(4.14) 



where A(9,9fi) is a non-negative function of its argu- 
ments. Here, Q has the same meaning that it had in Eq. 
(|3.ip - the rate at which heat is flowing into the system, 
in this case, from the reservoir - but now Q is a well de- 
fined function of 9, and Eq. (|4.14p is an equation, not an 
inequality. With this definition of Q, Eq. (|4.9p becomes 



9 Sn 



(4.15) 



V. SPECIFICS OF THE VACANCY MODEL 

Because we have an unambiguous definition of the 
total entropy, and because we know that the entropy 
of the vacancies must be included in it, we can write 
SneqiU, Vei,Nv) in the form 

SneqiU, Vel , N.,) = So{N,) + 5i(f/i) 

= SoiN,) + Si[U~ eo ~ C/ei(K/)] • (5.1) 

Equivalently, we can invert this relation and write it as 
an expression for the internal energy U 

UiSneq, Vel,N,) = C/o(7V„) + UiiSi) + UeliVel) 

= eoN, + Ui [Sneg - So{N,)] + UeliVel). (5.2) 

Here Uo{Njj) is the energy of the vacancies, eo is the 
formation energy of a vacancy, So{N.v) is the entropy of 
the vacancies, Uei{Vei) is the elastic energy, and Si and 
Ui are, respectively, the entropy and energy of all the 
other configurational, kinetic, and vibrational degrees of 
freedom in the system. The structure of these relations, 
i.e. the arguments oiUi and 5*1 in their second versions, 
describes the way the energy and entropy are shared be- 
tween the vacancies and the other degrees of freedom. 
Note that the total entropy and energy in Eqs. (|5.ip - (|5.2p 
are assumed to have very simple forms. For example, we 
have omitted a standard thermoelastic term proportional 
to SiVei in Eq. ([O]) . 

For specificity, we assume that the vacancies are very 
dilute, so that 

Sq{N,) = - TV, In [j^^ + N,, (5.3) 

where A'o is the number of sites at which vacancies might 
occur. Then, using Eq. (|5.2[) . we find that 



' dU 



[Ua{N,)~9So{N,)] 



eo -I- 6* In 



If we write 



— ) — 9 S„ eq = Cv f 



(5.4) 



(5.5) 



and interpret the extensive quantity Cv to be the heat 
capacity at constant volume, then Eq. (|4.15p becomes 



Cv9^ 
= Cv 



eo + P Wo + In 



No 



dC, 



(5.6) 



where the vacancy-related Gibbs free energy Gy is 

G,{9,p,N,)^eoN,-9SoiN,)+pvoN,. (5.7) 
The Clausius-Duhem inequality, Eq. (|4.12p . is 



Ny > 0. 



(5.8) 



Clearly, this term in Eq. (|5.6p is a non-negative rate 
of heat production associated with the relaxation of the 
internal variable Ny toward an equilibrium value. 



VI. EQUATION OF MOTION FOR iV„ 

At this point in the analysis, the standard procedure 
is to postulate a general form for an equation of motion 
for Ny, and to use the Clausius-Duhem inequality in Eq. 
(|5.8p to constrain the parameters that appear in it. In 
the present case, there is no reason why Ny should do 
anything more complicated than relax toward a stable 
equilibrium value. Therefore, for small departures from 
equilibrium, we write 



ToNy=fiNy) N^^9,p)-N, 



(6.1) 



where tq is a time scale, F is a positive, dimensionless 
rate factor that we anticipate will be a function of Ny 
(as well as 9 and p), and N^'^{9,p) is the equihbrium 
value of Ny at the given temperature and pressure. A 
convenient alternative form of this equation is 



..|^f(iV.) 



Ny 



f{Ny) ln(^ 



(6.2) 



Having no information about nonlinear corrections to Eq. 
(|6.ip , we can use the second expression on the right-hand 
side of Eq. (|6.2p just as well as the first, and will do so 
from here on. Other nonlinear equations of motion for 
Ny can easily be incorporated into this analysis when 
justified by some physical mechanism. 
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To satisfy the inequahty in Eq. (|5.8p . we require that 
both dGy/dNy, and the expression for Ny on the right- 
hand side of either Eq. (|6.1|) or Eq. (|6.2|) . vanish at the 
same point, i.e. at = N^'}{e,p). Thus N^i{e,p) is the 
solution of 



KdNy 



/ e,pM^=Nt,'' 



eo+pvo + e\n[^ ] =0, (6.3) 



and the equilibrium number of vacancies is proportional 
to a Boltzmann factor 



K9(0,p) =7Vo exp 



eo+pvo 



The inequality in Eq. (|5.8p is always satisfied so long as 
Ny is a monotonically decreasing function of Ny , which is 
required for dynamic stability, and is true for both Eqs. 
(lO) and (lO) . 

Retaining Q explicitly in this analysis has the added 
benefit of allowing us to deduce an expression for the 
rate factor r(A'^„). The equation of motion for Ny, as 
shown in Eq. (|6.ip . is a detailed-balance relation in which 
the vacancy creation rate is proportional to Ny'^{9,p). 
Therefore, according to Eq. (|6.4p . the creation rate auto- 
matically contains the appropriate Arrhenius activation 
factor, and F/tq can be interpreted as a dimensionless 
attempt frequency or, equivalently, a noise strength. 

Our experience with the STZ theory of plasticity leads 
us to write T as the sum of two terms: 



p{0) + riNy), 



(6.5) 



where p{0) is the strength of the noise generated solely 
by thermal fluctuations in the absence of mechanical de- 
formation, and r(7V„) is the noise strength associated 
with irreversible deformations, i.e. with nonzero Ny. For 
present purposes, we could simply set p{9) = 1 and let tq 
be temperature dependent; but we will need the explicit 
factor p{9) for discussing glassy systems in the following 
papers. 

An hypothesis (originally due to Pechenik [l^) that 
has worked well for the STZ theory is that F is propor- 
tional to the total rate per vacancy at which heat is gen- 
erated by the work done on the system. In the present 
case, this means that 



To 



■NyT{Ny)^9Sneq-Q = W, 



(6.6) 



where 9q is an energy, and W is the same non-negative 
dissipation rate that was defined in Eq. (|4.1ip . With 
this assumption, and with the second form of Ny given 
in Eq. (|6.2p . Eq. (|5.6p becomes 

- 9Ny\n^ (^^^ {p + T) + 9oTNy = 0. (6.7) 

Solving for p + F (a necessarily non-negative noise 
strength), we find 



T=p+T= 



l-i9/9o) In' {Ny/K^y 



(6.8) 



Thus, the mechanically generated noise enhances the rate 
factor, possibly quite substantially. The feature that F 
diverges when Ny is sufficiently far from its equilibrium 
value simply means that the system is dynamically driven 
away from such values of iVu, and that F remains positive 
at all times. 

Putting these pieces of the theory together, we have 



\n{Ny/N'^'i) 



Ny _ p{9) 

TQ l-i9/9o)ln{Ny/N^'')' 



(6.9) 



(6.4) Equation (|5.6p becomes 



Cv9 + 9 In 



Ny 



Ny 



A{9,9r){9r-9). (6.10) 



The combination of Eqs. (|6.9p and (|6.10p allows us to 
compute time dependent functions 9{t) and Ny{t) given 
any driving force p{t) and reservoir temperature 9ji{t). 
The simplest case is the limit in which the coupling to 
the reservoir is so strong that 9 — 9ii, and the heat 
capacity of the reservoir is so large that 9r remains a 
constant independent of how much heat is flowing to or 
from the system. This assumption, that the temperature 
is fixed by coupling to the reservoir, is implicit in most 
thermodynamic theories, but it is actually a bit subtle. 
The quantity Q on the right-hand side of Eq. (|6.10p is 
undefined in the limit A ^ oo, 9 ^ 9b,, which means 
that it can assume whatever value is needed in order to 
keep 9 = 9r = constant. With ^ = 0, Eq. ([eTO]) is just 
the Clausius-Duhem inequality again, now telling us that 
Q < 0, and that - as in the Kelvin-Planck statement of 
the second law - we are not allowed to convert heat 
directly into work. 

Alternatively, suppose that the process is adiabatic, i.e. 
A = Q. Then Eq. ()6.10p determines how the temperature 
rises as the external work is converted into internal heat. 
In the following paper, we will encounter an intermediate 
case in which A is small but nonzero. 



VII. CONCLUDING REMARKS 

Although the statistical analysis described here has 
been developed primarily for use in the effective- 
temperature theory of amorphous materials presented in 
the following paper [3| , the present results already point 
toward some general conclusions. 

We have shown that the Clausius-Duhem entropy- 
production inequality, when applied to the dynamics of 
internal degrees of freedom, can be derived directly from 
a statistical interpretation of the second law of thermo- 
dynamics - but only if the conditions listed following 
Eq. (|4.4p are satisfied. Perhaps the most important of 
these conditions is that the internal variables must be a 
small set of extensive quantities in order that the statisti- 
cal entropy be well defined for nonequilibrium situations. 
This condition, in turn, means that entropies as well as 
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energies associated with the internal variables must be 
included in any dynamical description of the system. 

The quantity appearing in the Clausius-Duhem in- 
equality in Eq. (|5.8p is interpreted by Lubliner as the 
"dissipation associated with the internal variables and 
their conjugate forces." Our analysis suggests a sharper 
and more physically intuitive interpretation - that the 
rate of energy dissipation W is the difference between the 
inelastic power —pVin = —p vo J^v and the rate of change 
of the free energy Uo{Ny) — 6So{Ny) that is stored in the 
internal degrees of freedom. 

One example of this difference occurs in Rice's 1971 pa- 
per [l3| , where he suggests that his internal variables rep- 
resent an extensive set of slips on slip planes distributed 
throughout a polycrystalline material. However, he does 
not calculate the number of ways in which the total slip 
can be realized as the sum of many individual slips, and 
therefore does not include the entropy associated with 
his internal "averaging variables" in dynamical formulas 
analogous to Eqs. (|5.7p and (|5.8p . 

These issues persist in the more recent literature and, 
in our opinion, are quite serious. For example, Anand 
and Su 18] implement something like Rice's picture of 
frictional slips on multiple slip planes by using a phe- 
nomenological, nonlinear, rate-dependent relation be- 
tween local flow and resolved stresses. Some memory 
of past deformation is carried by a "plastic volumetric 
strain" and by a related cohesion parameter that appears 
in the flow equation; but these are scalar quantities that 
cannot contain information about the directional history 
of shear flow. There is no dynamical yield stress as in 
STZ theory, nor - so far as we can tell - is there any way 
of using the theory to predict what happens when the 
loading stresses are removed or reversed, partly because 
the plastic volumetric strain and the cohesion parameter 
are scalars, but more importantly because neither are 
properly constituted internal state variables. The STZ 
theory has been developed explicitly to overcome such 



difficulties. It is the topic of the third paper in this series 
0. 

A related question, which is sometimes raised but not 
answered in the conventional literature, is what happens 
when the internal degrees of freedom do not relax in a 
simple manner. In the present case, our irreversible pro- 
cess can be described by a variational principle; that is, 
our single internal variable moves downhill in the one- 
dimensional free energy landscape defined by G.i,{Ny) in 
Eq. (jS.Zp . This picture is generalized in the conventional 
literature by assuming that a system with multiple in- 
ternal variables moves downhill in a multidimensional 
inelastic potential. The resulting fluxes obey what is 
called a "normality condition," or sometimes a "gener- 
alized normality condition" Q because they are assumed 
to be perpendicular to surfaces of constant (generalized) 
potential. 

We already know that the picture cannot be so sim- 
ple for the STZ theory, where increasing shear stress 
drives the system through an exchange of stability be- 
tween jammed and flowing steady states. This behav- 
ior is discussed in detail in the third paper in this se- 
ries [7]. More generally, we know that no such energy- 
minimization principles exist for many open situations, 
where the system is being persistently driven away from 
equilibrium, and where there are multiple, coupled, in- 
ternal state variables. It seems to us that it will be hard 
to predict a form for a generalized Clausius-Duhem in- 
equality without starting from a first-principles, fully sta- 
tistical and dynamical description of such systems. 
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